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Abstract 

We describe a linear homopolymer using a Grand Canonical en- 
semble formalism, a statistical representation that is very con- 
venient for formal manipulations. We investigate the properties 
of a system where only next neighbor interactions and an exter- 
nal, confining, field are present, and then show how a general 
pair interaction can be introduced perturbatively, making use of 
a Mayer expansion. Through a diagrammatic analysis, we shall 
show how constitutive equations derived for the polymeric sys- 
tem are equivalent to the Ornstein-Zernike and P.Y. equations 
for a simple fluid, and find the implications of such a mapping 
for the simple situation of Van der Waals mean field model for 
the fluid. 
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1 Introduction 

Molecules with linear chains as descriptives backbones can, in whole or in 
part, be constructed by temporally ordered sequences of additions. The ef- 
ficiency of such a process recommends chains as major actors in molecular 
biology, as well as in industrial chemistry. They may be examined at mul- 
tiple levels of resolution. In the study about to be reported, a very coarse 
description is used in which the units are taken as point entities with next 
neighbor strong binding forces and weak non-next neighbor interactions, all 
in the presence of a confining external field. Even more narrowly, we restrict 
our attention to homopolymers, but will indicate how questions concerning 
the statistics of heteropolymers can be addressed within this format. 

In further detail, we study a single homopolymer in classical thermal 
equilibrium. The homopolymer context suggests, in analogy with fluids, 
extension to a bath of monomers which can join or leave the chain, inducing 
fluctuations in the monomer number N and creating what might be termed 
a monomer grand ensemble. The corresponding formalism also follows the 
fluid model: suppose there are N' monomers in the full system, N of which 
are bound together to create a polymer, N' — N being free in the volume V 
envisioned. Then the canonical partition function of the complete system 
takes the form (one needs at least N — 1 to declare a polymer) 
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where Q^ is the monomer canonical partition function in its center-of-mass 
coordinate system, Q&> that of the iV-monomer polymer in the external 
potential field u. The "grand" partition function of the monomer collection 
defining the polymer will then be taken as 

SW(C)= hm — N , Q T N % (2) 

or via the Stirling approximation, 
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identifying the fugacity as 

£( r ) = Coe-^W = e^ (r) . (4) 

The expression Q is very convenient for formal manipulations, but a certain 
amount of caution must be exercised. In particular, number fluctuations 
can be very large, as we will see, so that "canonical" - i.e. fixed monomer 
number - information, requires care. 

In this paper, our aim is to show how (jSJ) allows for a familiar type of 
diagrammatic representation, which not only gives rise to a systematic per- 
turbation expansion in the strength of the non-next neighbor forces, but 
also - under special but physically reasonable restrictions - maps onto a 
suitable classical fluid, with its array of well-tested computational recipes. 
With this, we can routinely examine the effect of structured fields on the 
polymer density profile, which is our implicit objective. 

2 The Reference System p 

Now let us introduce explicit notation following [2j. The model homopoly- 
mer that we will deal with has as monomer units point "particles" located 
at {vi,i = 1, . . . , iV}, which may be imagined as spatial locations, but can 
include type, internal state, etc. There is an external potential field u(r) 
appearing only in the combination £(r) = Co e~^ r \ where (3 is the recipro- 
cal temperature, and a next neighbor potential <ft(r,r'), typically <p(r — r'), 
with associated Boltzmann factor w(r,r') = e~ l3 ^ r,r ' as well as a general 
pair interaction (including next neighbors) v(r, r'), giving rise to the Mayer 
function f(r, r') = e~^ r,r •* — 1. Note that we can expect that w — > 0, / —> 
as \r — r'\ — > oo. When r represents real space location alone, we have the 
format of freely jointed units, but that is not intrinsic, since the units can 
involve more than one monomer (see e.g. [3]). 

We will choose, as generating function for the thermodynamics, the grand 
canonical potential 

fi[C] = -(l//?)lnH (p) [C], (5) 

so that , e.g. the resulting monomer density profile will be given by 

n(r) = -C(r) w 

= + ^cw (6) 



To start with, consider the "reference system" in which the general pair 
interactions have not yet been included. Since the monomers can be labeled 
ordinally along the chain, and they are not equivalent with respect to system 
energy, there is no statistical weight factor for an iV-monomer configuration. 
Hence (JHJ) can be written as 
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where w is regarded as a (continuous as needed) matrix, C a diagonal matrix, 
/ (I(r,r') = S(r,r')) the identity, and 1 the vector of all l's. The validity 
of (jZJ) is subject to convergence. In terms of the molecular weight control 
parameter ( , it is apparent that (JJJ) will converge only until Co max = 
^-/Imax (we _/3u ) - the value at which the minimum eigenvalue of w ( reaches 
1. The behavior of the system as Co "~ > Co max must be factored out. 
Combining (JHJ) and ((7|), we have, using the standard (A^ 1 )' = —A^^^A'A -1 
for any derivative (') on any matrix inverse, the mean density profile 

nM (1|C (/ - C^)' 1 k)C(r)(r| (7 - ^Q' 1 |1) ( „s 
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In particular, the mean number of monomers, iV = J n(r)dr, becomes 

w= (nai-vcr 2 \i) (9) 
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dominated by iV = Co max /(Co max - Co) as Co -> (max- Thus, large mean 
molecular weight occurs very close to the "resonance" at C = Co max- The 
effect is brought out forcefully by computing the standard deviation of N 
in a similar fashion: 

6N- 1/2 
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evidence of something like an exponential distribution. 



Figure 1: Example of diagram in expansion (|llj) . 



It is also useful to look ahead and see how the qualitative behavior will 
be modified in the presence of non-next neighbor interactions, which of 
course contribute N 2 terms to the configurational energy. If all N 2 were 
positive and of similar magnitude (e.g. long-range) they would append a 
convergence factor of the form e~ aN to the series (J7J). This strategy is being 
investigated. However the non-next neighbor forces of physical interest are 
not necessarily long-range, and have a substantial attractive component. 
If the monomers are prevented from clustering, their forces would again 
contribute to the order of N, and would simply modify the reference system 
pathology. Short range forces would reinforce this prohibition, but in the 
model we will attend to, it is mainly the profile shaping due to the confined 
external potential that is relied on. 



3 Diagrammatic Analysis 

We proceed to the fully interacting homopolymer ensemble. In terms of 
e -/3v{r,r ) _ y _|_ f( r ^ r ') 5 the N th term in (J7|) becomes 
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(11) 
where II' is the product over ^N(N — 1) distinct unordered pairs; thus a 

typical term of the full expansion, in terms of multinomials on f, would be 

diagrammatically like the one shown in FigfTTl with multiplicative Boltz- 



mann factor represented by: 

• = C(r) external field 



w(r,r') next neighbor weight 
/ (r, r') general pair Mayer factor 



(12) 



All field points r* are integrated over. It will be convenient to introduce the 
matrix: 0*wv^ww*wO denoted by A, for N > 2, in which the endpoints 
have weight 1 and are not integrated over, but the sum of all diagrams over 
all N > 2 is included. Hence 

2 (P) [C] = (1|C) + (C|A|C) (13) 

where the vector \() = (\l). 

Now various diagram reductions can be carried out. The basic one stems 
from the observation that some of the vertices are articulation points in the 
sense that their removal disconnects the diagram into a left and right half, 
e.g. r 6 and r 9 in FigHJ 

The end points are never given the status of articulation points. If the 
matrix belonging to the sum of all connected diagrams - those with no 
articulation points, those with endpoints (O) unintegrated, is denoted by 
C, then dearly A = C + CCC + CCCCC + ..., or: 

A = (I-C(y l C } (14) 

and (JT3j) can be written as 

H^[C] = (1|C((/-CC)" 1 |1)- (15) 

Of course, in the absence of the general pair interaction, C = w, and 1)15)1 
reduces to (J7J). 

But C can also be built up from A. To do this, observe that C consists 
of two types of diagrams, C%, in which the endpoints have a direct f-link to 
each other, and Cq in which they don't. See FigElfor examples. Since the 
end-linked diagrams Ai of A are precisely the end-linked diagrams of C\ of 
C, and since C\ — f * (A — Ai) we have 

C 1 = f*(A-C 1 ), (16) 
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Figure 2: Examples of graphs included in A,Ci, and Co- The shaded gray areas 
indicate the presence of any possible subgraph. 



where A * B denotes element by element multiplication of the matrices A 
and B. On the other hand, Cq is composed of w, together with all end- 
unlinked diagrams that have no articulation points. The latter are built by 
inserting A's between the termini of any internal f links, and so 



Cq = W + O v ^^^# WAWv # w ^^O + 0*W*M%H*W+^M*MMM%fJAM+,0 + . . . (17) 

But (jTrT|) can be written as 

e*Ci = /*A, (18) 

where e = 1 + /, leading us to the second relation between A and C: 



e*C = e*w-\- O^ww'O + e * o^^ w #^^^ w # ww ^ v O + • • • (19) 
to accompany (|14|L now written as 

A _(7 = CCA (20) 



4 Mapping onto Fluids 

The practical theory of classical simple fluids in thermal equilibrium has 
had a long period of development, entailing a sequence of more and more 
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sophisticated approaches jlj. For uniform simple fluids not at thermody- 
namic singularities and not with long range forces, little remains to be done. 
And the non-uniform situation is increasingly under control. Let us indicate 
that we are speaking of a fluid - initially related to the homopolymer under 
study only by having a common pair interaction, Vl(t, r') = v(r,r') - by 
appending a subscript (L). For present purposes, it's sufficient to observe 
that almost all of the older approximations posit an algebraic relationship 
(or functional relation in more recent versions) between the dimensionless 
pair correlation function 

Mr>0 = -25^-1 (2i) 

n L (r)n L (r') 

and the direct correlation function CL{r,r'), directly related to the external 
potential - density linear response. 

In classical thermal equilibrium, hi and cl satisfy identically the Ornstein- 
Zernike equation [5], written in continuous matrix form as 

h L -c L = c L n L h L , (22) 

where n^ is the particle density as diagonal operator :(r\ni\r') = riL(r)S(r — 
r'). We then need a second relation or closure to obtain both hi and cl- 
One of the most effective reasonable approximations for short-range forces 
is the P.Y. relation [6], which in terms of Boltzmann factor e and Mayer 
function / = e — 1 of the pair interaction v reads 

e*c L = f*h L + f, (23) 

where again the (*) denotes element by element multiplication. 

Comparison between the pair (|22I23|) and the pair (|2UI19|) . is suggestive. 
If we terminate (|19j) at first order: 

e*C = f * A + e*w, (24) 

becomes identical with (}2"2"j) under the assumption 

C = Kc L} k = Kh L} ( = n L /K (25) 

if the condition, for any fixed constant K, 

Kf = e + w (26) 




Figure 3: Example of eq.(|27|) 



is satisfied - and (J2*U|) and (|2*3*j) are identical as well. Expressed in terms of 
the general pair interaction v and the next neighbor potential 0, eq. (ffilj) 
says that 



+ 



K 



is satisfied. For example, it is sufficient to choose v as, qualitatively, a 
reduced version of the highly attractive localized <fi. Note that according 
to our correspondence, m = K( may appear as a quite high density in its 
confined domain. Clearly, (1//3) InK + (f) > is required by (|27jl . 

Once we have used the P.Y. approximation to produce a fluid correspond- 
ing to the homopolymer under consideration, we can, in a logically tenuous 
fashion, analyze the fluid by any other convenient approximation method. 
Whatever approximation is used, the derived end result is the expression 
of the polymeric monomer density resulting from the imposed contain- 
ment potential represented by C( r )- It's only necessary to translate this 
into equivalent fluid language. From JBJ), (|T3|l and (|25|). we have at once 
(1|C + C AC|1> = i(l\n L + n L h L n L \l), or 



•s(p) 



s(p) 



(l\S L \l) 
K 



(28) 



where Si(r,r') = 715(7", r') — ni(r)riL(r') + riL{r)S(r — r') is the complete 
Ursell function or generalized structure factor of the fluid, also appearing 
as the density- density correlation function. 



5 Discussion 

The picture of an equivalent non-chain fluid is a bit deceptive. At first 
glance, it resembles the familiar statistical model of biopolymers, itself an 
offshot of the Lifschitz homopolymer condensation transition model jjj. But 
the difference is seen most vividly in a situation in which there is an un- 
constrained volume - constant external potential - bounded by a hard wall 
container. The "dual" fluid to the polymer would then have a constant 
density, terminated by zero density at the boundary, and an internal po- 
tential structured to produce the required uniform fluid. Some meaningful 
conclusions can be arrived at in this situation without an in-depth analysis. 
Suppose we restrict attention to observations at the correlation length 
scale of resolution. On this scale, one knows that, clearly, 

S L (r,r') = S L (r-r';n L {r)) 

= 8{r-r') [ S L (R;n L (r))dR (29) 



where Sl(R',til) is the density-density correlation for a uniform system of 
density til- Since 

S L (R;n L )dR=^, (30) 



this implies that 



and consequently that 




S L (r.r')dnlr' = / J^W, (31) 
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n{r) = C(r)^jlnE«[C] 

5 5n L (r) f 6n L (r') , 
= n L [r)- — — — — — / / — — dr, 32 

5n L (r) 5f3/j, L {r) J 6(3 n(f) 

all of which is at ni(r) = KQq. This would seem to indicate that the floppy 
chain, on averaging over all configurations, would simply fill available space 
uniformly (since n^(r) has no spatial dependence). But this is only true if 
the correlation length is much smaller than the diameter of the confining 
volume. For example, in the primitive Van der Waals mean field model |H] , 
in which 

lnn L (r) = [3fi L (r) + i f n L {r)n L {r') ( e -/M»v / >- 1 ) dr >, (33) 



it is easy to see that (J52J1 for a uniform fluid would read 

1 fl-2n 2 T B\ . . 

where B is the second virial coefficient, negative for attractive v, and V 
the confinement volume. Then, the needed divergence of n V for large V 
at "resonant" n, would imply a divergent correlation length, and so in an 
exact solution, the density would rise only gradually from its surface to 
bulk value, compatible with a confined polymer transported throughout 
the volume. 

One can of course also go beyond the equivalent fluid picture, but certainly 
the domain untouched in our treatment is that of heteropolymers. Statisti- 
cal models are readily constructed by including species dependences in the 
one and two-body potentials, but the statistical distribution is broad unless 
the units are at least taken as short subsequences of monomers, thereby 
introducing multiple internal degrees of freedom. This strategy is being 
investigated, and in fact was implemented at a very crude level some time 
ago in the context of including bond angle and dihedral angle restrictions 
to accord with reality even in heteropolymers [3]. 

6 Acknowlegdments 

S.P. would like to thank PCT-ESPCI and Volkswagenstiftung for current 
support, and the New York University Department of Physics for the sup- 
port during her Ph.D. when this work was mainly conducted. The work of 
J.P is partially sponsored by DOE, grant DE-FG02-02ER1592. 

References 

[1] See also H.L. Frisch and J.K. Percus, J. Phys. Chem B, 105, 11834 
(2001) 

[2] J.K. Percus, J. Stat. Phys, 106, 357 (2002) 

[3] K.K Muller-Nedebock, H.L. Frisch, J.K. Percus Phys. Rev. E, 67, 
011801 (2003) 

[4] Hansen J. P. and McDonald I. R., Theory of Simple Liquids, Academic 
Press, London (1986) 

10 



[5] Ornstein and Zernike Proc. Acad. Sci. Amsterdam, 17, 793 (1914) 
[6] J. K. Percus and G. J. Yevick, Phys. Rev. 110, 1 (1957) 
[7] I. M. Lifschitz, Sov. Phys. JETP, 55, 2408 (1968) 
[8] J. D. Van der Waals, Z. Phys. Chem, 13, 657 (1894) 



11 



